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Abstract: Thermostats models in space dimension d = 
1,2,3 for nonequilibrium statistical mechanics are consid- 
ered and it is shown that, in the thermodynamic limit, the 
motions of frictionless thermostats and isokinetic ther- 
mostats coincide. 



I. THERMOSTATS 

A test system of particles in a container fio a-nd v sys- 
tems of particles in containers fli, . . . ,^1^, interact and 
define a model of a system in interaction with v ther- 
mostats if the particles in fti, . . . can be considered 
at fixed temperatures Ti , . . . , . 

A representation of the system is in Fig.l: 





= (Xq, Xq, Xi, Xi, . . . , Xjy, Xi,) 

Fig.l: The 1 + v finite boxes Qj H A, j = 0, . . . , v, are marked 
Co , Ci , . . . , Ci/ and contain No, Ni, . . . , N^, particles out of the 
infinitely many particles with positions and velocities denoted 
Xo, Xi, . . . , X^, and Xo, Xi, . . . , X^/, respectively, contained 
in ^Ij, j > and considered in the d — 1,2 cases. The second 
figure illustrates the spectal geometry considered for d — 3 (as 
well as for d = 1, 2): here two thermostats, symbolized by the 
shaded regions, fii , ^2 occupy half-spaces adjacent to f2o . 

A formal description of the model can be found in [1,2]. 
The temperatures in the thermostats will be identified 
with their average kinetic energy per particle. 

Considering external thermostats as correctly repre- 
senting the physics of the interaction of a system in con- 
tact with external reservoirs has been introduced in [3]. 
Their analysis was founded on the grounds of 

(1) identity, in the thermodynamic limit, of the evolution 
with and without thermostats 

(2) identity of the phase space contraction of the ther- 
mostatted systems with the physical entropy production 
(up to a time derivative). 

Here we follow their strategy. To implement the phys- 
ical requirement that the thermostats have well defined 
temperatures and densities the initial data will be imag- 
ined to be randomly chosen with a suitable distribution. 
If A is a ball centered at the origin the space 'H(A) 
will be the space of the finite configurations (X, X) with 

X c Uj>of^j n A. 



Referring to Fig.l denote by Kj{ii.j) — ^X^ the kinetic 
energy of the particles in the j-th container (assuming the 
particles mass = 1), by Uj{'K.j) their potential energy 
(more explicitly expressed in Eq.(1.3). 

Initial data: The probability distribution for the random 
choice of initial data will be, if dx 11^=0 '^'^'tf^' ' 
distribution on 'H(A) 



fj.A{dx) 



^n/(^M(^)-4iv.)rfx (1.1) 



wtth Ho{x) = T:,=oPj{Kj{±,) + C/,(X,)), and N, = 
(5j I A n Vlj I with Pj T^jr > O7 '^i > 0, j > ; the values 

Pa = 'k~To > '^0 > will also be fixed (but bear no par- 
ticular physical meaning because the test system is kept 
finite). 

Here 6 = {Sq, Si,... 6^) and T (Tq, ^i, ■ • ■ , T^) are 
fixed densities and temperatures, A is a ball centered at 
the origin and Z is the normalization. 

With the above choice of initial data the physical re- 
quirement that the thermostats are in a configuration 
with densities and temperatures assigned is realized at 
time 0. But the probability distribution is not invariant 
under time evolution: because the temperatures f3j are 
different and also because Hq does not contain the in- 
teraction between the system and the thermostats and 
because on the system will be imgained to act noncon- 
servative forces. 

The equations of motion (see Fig.l) in A = A„, with A„ 
being the ball of radius 2'^r^ with a length unit that 
can be chosen to be the interaction range (see below), 
will be different in the frictionless thermostats and in 
the isokinetic thermostats. They will differ by the value 
of a parameter a = 0, 1: 



mXo. = -a,f/o(Xo) - E,>o drUoASia, X^) + *,(Xo) 
mXji = -dtUjpi-j) - 9it/oj(Xo, Xj) - aajXjj(1.2) 

where the first label, j = or j = 1, . . . ,v respectively, 
refers to the test system or to a thermostat, while the 
second indicates the components of the coordinates of 
the points located in the corresponding container and 
initially in the regularization box A„ (hence the labels i 
in the subscripts {j,i) have Nj d values). Furthermore: 

(1) the #(Xo) are, positional, nonconservative, smooth 
stirring forces, possibly vanishing; 

(2) other forces are conservative and generated by a pair 
potential with range r^p, which couples all pairs in the 
same containers and all pairs of particles one of which is 
located in fig and the other in Q,j [i.e. there is no direct 
interaction between the different thermostats). 

(3) particles are repelled by the boundaries d^j by a con- 
servative force of potential energy ijj, of range <C r^. 
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diverging at the walls. The potential energies will be 
UjCKj), j > 0, and C/o.j(Xo, Xj): respectively denoting 
the internal energies of the various systems and the po- 
tential energy of interaction between the system and the 
thermostats. 



C/o,j(Xo,Xj) = Egexo.g'ex, - 1') (1-3) 

(4) in the case a = particles will be allowed to exit 
the regions Qj fl A„ through the boundary of A„ while 
if o = 1 they will be constrained to remain inside the 
container ilj O A„ by an clastic reflection on (9A„: this 
choice of the boundary conditions is imposed mainly to 
make possible references to [2] . 

(5) the aj are determined so that the solution to the 
equations with a = 1 keeps exactly the same kinetic en- 
ergy it has at time 0, i.e. ^NjkBTj, for j > 0, see 
Eq.(l.l); this means 



(1.4) 



where Q 



def 



X, • 9x, 



t/oj(Xo,Xj), is the heat ceded 



by the system to the j-th thermostat. 

(6) The potentials ip, tJj have been chosen j~independent 
for simplicity. The pair potential (p will be supposed 
smooth, > 0, with finite range r^p > 0, setting (p{0) 

def 

= (fo > 0; the wall potential will also be supposed 
smooth at a distance r > from the walls, > 0, with 
range r^; , and diverging proportionally to r~" as r — )• 0, 

for some a > 0. 

Hypothesis: The initial state distributions /xa„ are as- 
sumed to satisfy a large deviations property for the total 
potential energy f/j,A„(a;); in the sense that for C > 
there is c> such that: 



MA, 



l^^j n A„| 



> 



}) 



< e 



-c2 



l(l-2E)d 



(1.5) 



for suitable Uj G M, e G [0, 1). This is a "no phase transi- 
tion" assumption ( satisfied in the cluster expansion region 
of the parameters T, A as discussed in [4, 5]). Notice that 



the volume 



vtj n A„ 



is 0(2"''). 



The analysis will be restricted to the geometries in 
Fig.l if d = 1, 2 and only in the second of Fig.l if d = 3. 



II. TIME EVOLUTION 

Infinite systems are idealizations, not uncommon in 
statistical mechanics, that must be considered as limit- 
ing cases of large, yet finite, systems. We therefore call 



regularized equations of motion the Eq.(1.2) and denote 
S'j"'°^x, or also x^^'°'\t), their solutions with data ran- 
domly chosen with the distribution in Eq.(1.2). 

Solutions exist (for all data in 7^ (A), see the initial 
data hypothesis in Sec. I) by the standard existence and 
uniqueness theorems for ordinary differential equations 
in the case a = 0. 

The case a = 1 involves clastic reflections on dK and a 
proof of existence of the evolution requires (by adapting 
the analysis in [6], see [2, Appedix I]) obtaining exixtence 
and uniqueness for all initial data outside a set of fi\^ - 
probability 0. 

We shall discuss results that hold uniformly in the size 
n provided it is large enough. The results can be conve- 
niently formulated in terms of quantities defined below. 

Let la;. _ ^'1 %f Mi + taU, and 



+ k E,y/» - Ij) + + Vq) (2.1) 

Ni(x) number of particles within r^ of qi, 

def 



Let log+^; = max{l,log2 |^;|}, g,:{z) = (log+2;)^, p„(g) 

def 

= distance of q from [Jj{dClj U 9A„) and call Xe = 
{x \£{x) < E} with 



c/ \ def 

t[x) = sup 



sup 

? R>gi/d(x/r^) 



W{x;^,R) 



(2.2) 



The set Xe has /Ua„ -probability approaching 1 as ^ 00 
uniformly in n (see, for instance, [2, Eq.(9.2)] choosing 
c = E and c' = l,7(c) = c — Cq): 



I^aAXe) >l-Ce 



-cE 



(2.3) 



for suitably chosen C, c > 0. 

Let d < 3, then given an initial datum x the motion 
exists for /iA„ almost all x e Hi/a and \/t > 0. 
Fixed arbitrarily an observation time 9 < 00 the main 
result will be 

Theorem 1: // x G 7^(A„) n Xe and < t < B there 
are c = C{E,e) < 00, c' = C'{E,Q) > 0, 1 > 7 > 1/2 
such that, for particles which at time are in A^, i.e. 
qi{0) € Afc. Consider, for n> k, the events 



(1) 
(2) 
(3) 

(4) |(5("'°)x),-(5, 



q'r'°Ht)\<cv,k^, 



Pn{q^'''\t))>c'r^k- 



+i'°)™),| < e-'2"^^ n>k, 



(2.4) 



(5) |(5t^"'°'a;)i - {S^'''^^x)i\ < e-^'(i°s")^ k < (logn)T. 
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The events (1-4) are realized for all x e Xe, while 
the event (5) is realized with -Probability iTn and 
|/^A„('^f;)-vr„| <ce-='(i°sn)=\ 

This means that, if A„ is large, motion of the par- 
ticles close to the test system is largely independent on 

the rcgularization size n. And thermostatted motion and 
frictionless motions, near the test system, arc also very 
close. 

The uniformity in n of the constants c, c' is the really 
interesting part of the statement. 

Items (1-4) arc proved in [2, theorem 5, sec. 6] for d = 
1,2,3 (and for d = 1,2, only, in [1, theorem 7] via a 
different method). Therefore the novelty in the present 
paper will be the proof of item (5): it will heavily rely on 
[2, Sec.VII]. 

For later reference it is convenient to introduce a few 
more notations. 

Let the cube with side r^ centered at a point ^ in 
the lattice Ti^Z**. Let A be a subset of R** and, using the 
definitions in Eq.(2.1), 



def 



def 



NAixy^ EgeA^, l^A"=max,eA^ (2.5) 



II II 



max(iV<7^(a;),ec^ (a)) 

Si/2(l/i-v>) ' \ ■ ) 



see Eq.(1.5) and Eq.(2.6) for notations. Split Xe = An U 
Bn with 



Bn''={xeXE:Tn{x)<e}: 



(3.4) 



Theorem 2: There are positive constants C,C',c de- 
pending only on E such that for all n large enough: 

(1) if t < Tn{x), Sf^'^^x and S^^'^^x are close in the 
sense that for qi{0) G ^{\ogn)i 

\q^''\t) - q^'°\t)\ < Cr^e-('°s")^^ 
\Qi"'^\t) - Q^'°\t)\ < Cvt e-a°s")"^ (3.5) 

Furthermore with the notations in Eq.(3.1), (3.2), for n 
large enough and for all t < Tn{x): 

Na, (5t^"''^a;) < C*, ^A. (^l"''^a:) < C* + 1 (3.6) 

(2) the set Bn has iJ,An~probability bounded by 



(3.7) 



where ec^ix) ^eQ(a;), ec{x) = maxggc(^9^ + V'(a))- 



III. THERMOSTATTED EVOLUTIONS 

Thermostatted evolution can now be studied by com- 
parison with the frictionless one. 

The problem will be studied by restricting attention to 

a suitable subset of the set Xe =■ \x\E{x) < 

Consider the bands of points ^ at distance po(C) within 
r^ and 2r^ from the boundary dflo of SIq- 



A* =^ {q : Po{q) < r^}, A,, 1?/ 1^ . ^^(g) < (3.1) 

By items 1,2,3 theorem 1 in Sec.II (which can be taken 

for granted by the remark following theorem 1) there is 
C* > (depending on E) so that, for all x G Xe and 
with the notations Eq.(2.6), for n large enough: 



max max| A^A., (Sf 

i<0 



^"'°^a;), FA..(5^°^a;)}<a (3.2) 



Fixed 7 once and for all, arbitrarily with i < 7 < 1, let 
S' be the set where are realized the events delimiting the 
stopping time Tn{x) 



T„{x) { maxi : t < 9, Vt < i. 



(3.3) 



|n,nA„| 



— 



<||,||5i"'^MU<(logn)^}. 



Remarks: (1) The Eq.(3.5) together with the first four 
items of theorem 1 imply Eq.(3.6). 

(2) A sequence of initial data {a;} = {xi, X2, ■ ■ ■) sampled 
randomly and independently with the distributions /Ua„ , 

i.e. with the distribution fio(d{x})'^= Jl^i MA„ ('i2;„), 
will consist of configurations Xn ^ Bn for all n large 
enough, by Borel-Cantelli's lemma and Eq.(3.7), because 
7 > 1/2. 

(3) To prove item (1) and Eq.(3.6) we shall compare the 
evolutions x^"'^'>{t) with x^"'^\t), at same initial datum 
X G Xe and t < Tn{x), the latter being the stopping time 
defined in Eq.(3.3). 

Two preliminary results are necessary, namely that 
there is C > so that for all n large enough the fol- 
lowing holds. 

Lemma 1: Let x € Xe, t < Tn{x), see Eq.(3.3), 

k > (logn)''', then 



ii{t)\ <Cvx{k\ognf 



9l"'''(t)| <r^(2*= + C(A;logn)T). 



(3.8) 



for qi{0) e Ah and t <Q. 



This lemma is needed because, otherwise, the positions 
and speed at time t cannot be controlled in terms of the 
norms ||x||„ at time 0: since the particles move they must 
be followed (a "Lagrangian" viewpoint). 

A corollary of the above will be: 
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Lemma 2: Let M and p he the maximal number of par- 
ticles which at any given time < T„(.'r) interact with a 
particle qi initially in A^^i and, respectively, the mini- 
mal distance of a particle from the walls. Then 

M < C [klognf'^ , /)> c(fclogn)-2(<^T+^)/" (3.9) 

for all integers k e ((logn)''', 2(logn)'>'). 

Rem,arks: (1) a is the power which controls the diver- 
gence rate of the wall potentials. 
(2) The proof of the lemmas is in [2, Appendix L]. 

Proof (of theorem 2): The proof of the key "entropy 
bound" Eq.(3.7) is similar to the proof of the correspond- 
ing statement in [2] and is reproduced in the Appendix 
A, because it requires some changes with respect to the 
analysis in [2] where the stopping time definition was 
based on the kinetic energy rather then on the potential 
energy as in Eq.(3.3). In this appendix we make essential 
use of the hypothesis at the end of Sec.l. 

We shall now bound 6i{t,n)'^= |gf ''^(i) - 
again following [2]. Let fi be the acceleration of the 

particle i, with qf^'"""^ {t}) G A^+i, t < r„(x), and k'^= 
(log_|_ n)'', due to the other particles and to the walls. 

If Qi = qi{0) G Ak, \fi\ < C (fclogn)"', v'=^{2dj + 
1) (1 + ^) by Eq.(3.9). Consider the two evolutions, for 
a = 0, 1 respectively, 

where the label j on the coordinates (indicating the con- 
tainer) is omitted and fi is the force acting on the se- 
lected particle divided by its mass (for j = it includes 
the stirring force). 

Subtracting the Eq.(3.10) for a = and a = 1 it follows 
that for any qi G Afc (possibly close to the origin hence 
very far from the boundary of A), if n is large, because 
k = (log+ n)T) 

Si{t,n) <C{k\ogn)'^'2-"^^ (3.11) 
+e/ol/i((?("'iHT))-.A(g("'°H^))|rfr. 

provided \ J^^'^ ajds] is bounded proportionally to 2""^'' 
for [ii,i2] c'[0,T„(a;)]. 

Prom the definition Eq.(1.4) of aj the bound of the 
integral \ J*^ C(jds\ with [ti,t2] C [0, T„(x)] is split into 
a bound on Jg^'^'^^Sj^-^^dt and a bound on the ratio 

Pj dNj 

The first bound can be derived from the inequality 
||S'i"'^^x||„ < (logn)T [i.e. for t < T„(a;)) and it is < C 



(fclogn)'''2~"'^. This requires using lemmata 1,2 because 

|jS't.T||,i only gives information about the particle that at 
time are close to the test system, while the i-particle 
might be close to Cq at time t but not at time 0. 

The second bound also follows from the definition of 

the stopping time, which implies the validity of the in- 

equality "'^'^'''^"'"y < C2--'^. This also 

makes use of lemmata 1,2 (for the same reason as above). 

Let £ be a non-negative integer, k^ such that 

2^' =2'' + (.C{k\ogny (3.12) 

and Ukt{t,n) the max of Si{t,n) over \qi\ < 2'^^. 

The difference in the accelerations is bounded, by 
lemma 2, by the maximum number (fc log n)'''^ of the par- 
ticles which can interact with qi(t,n) times max|5^y)|, 
plus a term proportional to (k\ogn)'-^'^''^^^'-^^^^"'> due 
to the walls potential. Then by Eq.(3.11) and writing 

r?"1i/(2rf7+l)(l + ^), 

< C(fclogn)'''2-"='* 

+C(fclogn)"'7o'^i^f. (3.13) 

for £<t = 2^/{{k log nyC), the latter being the largest 
(. such that 2'='^ < 2^+^. By Eq.(3.13) 

Uk{t,n) < eC'C' i°s")''"c(A:logn)'''2-"^^ 

^ (ccfciognr-r c{k\ogn)^. (3.14) 

Thus for n large enough Uk{t,n) is bounded by the 
r.h.s. of the first of Eq.(3.5); analogous argument shows 
that also the velocity differences are bounded as in 
Eq.(3.5) which is thus proved for all t < Tn{x). 

It follows that, for n > e o and i fixed, given 
<Z,(0) with \qm\lr^ < 2>'o it is \ql"''\t) - qt"\t)\ < 
U(\ogn)i(,t,n) < C'e~'=('°s"^^, i.e. for n large is 
closer than r^ to qi"''^\t). 

Remarking that we know "everything" about the 
Hamiltonian motion we can use such knowledge by ap- 
plying Eq.(3.5) to particles which are initially within a 
distance r;^2'^° of the origin, with ko fixed arbitrarily, for 
all large n. 

Therefore the number of particles in ^{"'^^(f) which are 

in A* is smaller than the number of particles of q^"'^^ {t) 
in A** which is bounded by C*. 

An analogous argument for the velocities allows to 
complete the proof of Eq.(3.6), given the closeness of the 
positions and speeds of the motions with a = 0, 1, see 
Eq.(3.5). Hence T„(x) = © unless x G Bn and the proof 
of theorem 3, hence of theorem 1, is complete. 
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IV. APPENDIX A: ENTROPY 

Entropy production per unit time, in a configura- 
tion X, is naturally defined in terms of Qj = — Xj • 

i^Xj t^Oj (Xq, Xj), interpreted as the heat ceded by the 
system to the thermostats and it is given by: 

ao{x) = Y,p,Q,{x) (4.1) 
j>o 

In the isokinetic thermostat model, a = 1, and if the 

volumes in phase space are measured by the distribution 
/Ua„j a direct computation shows that this quantity dif- 
fers from the contraction rate of the phase space vol- 
ume by MQo + $(^o) • ^o) = + Uo) (Qo being 

— Xq • 3xoC^,j(Xo, Xj)) and by a further "small correc- 
tion" see also [7]. It is: 

a{x)=ao{x)-^P,^^^^+Po{Ko + ifo)- (4.2) 
j>0 J 

The Eq.(3.6) yields a bound |ao(a;)| < C, ^ < 

(72"""^ < C and also the last term in Eq.(4.2) gives a 
similar bound (because the number of particles in Hq is 
bounded in terms of and Kq+Uo equals (3o+^(Xo)-Xo 
is bounded by Eq.(3.6) and, via the comparison with the 
frictionless motion in Eq.(3.5) and the first four items in 
theorem 1). 

Furthermore by the definition of the stopping time it 
is J*^ w-j'^t< C2-"^'^: so that there is C and 

r \u{S'f-''^^ x)\dt < C , yti<t2<T„{x) (4.3) 

Writing fcj for the smallest integer > {log n)^gj{^/r,j)) 
(here is chosen instead of the natural 51/2 in order to 
simplify the formulae: recall that by defiinition 7 > 1/2), 
then codimension 1 surface containing S' splits into an 
union over ^ e A„ n r^Z"^ of the union of 5| U <S| U 5^, 
where 

4= {yeJ:' ■.\ync^\ = k^,\yndC^\ = i} 

5| = {y e E' : y n Q 9 {q, q), e{q, q) = 

4= {y€^':!^^-u,=±^} (4.4) 

if £| =^ ((logn)'>'5^(C/r^))^ and K Qj n A„. 



By the assumption Eq.(1.5) it is not restrictive to sup- 
pose that att = it is | ~ < 52^}- 

Abridge Aj^n — Uj and consider first the 

case of . Let D C be the set of the x which satisfy 
Aj n{x) = (^o2~"^'' for a given j > while Aji nU{x) < 
(^o2-"^'^ for j' > 0, 3' ^ J. 

Then the probability ij.\^{D) can be bounded, see 
Eq.(3.2), by 

e^' I MA„ {dx) e S{Aj,n - G) \G\ x{N < p2"'') (4.5) 

where G is any value between and and G{x) 

is the derivative Aj^a i.e., from the equations of motion, 

G{x) (E(9 - ^')9,'P{<1 - 9') + E ^^-^ (4.6) 

" 9.9' 9 

and e*^ takes into account the entropy estimate i.e. the 
bound Eq.(4.3) of the non-invariance of /iA„- 

Averaging the bound in Eq.(4.6) over G between \ ^p-\ 
and our bound is 

I I 

e^'l^ / MA„(da;)|G(a;)|xx (4.7) 

where Xx is the characteristic function of the set where 
€ |v^' Jy^]' Schwartz's inequality this is 

<e^'^(MA„(G2)^/VA„(x.)^/^) (4.8) 

By superstability and the integral can be bounded above 
by C while the second integral estimates the square 
root of a probability of a lage deviation and is there- 
fore bounded by by e^*^^ . Hence the contribution 
of C3 to the boimd on /iA„ {B) is (amply) bounded by the 
r.h.s. of Eq.(3.7). 

Similarly, the surface areas /iA„.S'('5|) and /iA„,S'('5|) 
induced by /xa„ on 5^,5^ induced by no are bounded by 

MA„,S'(5|) < Cv^e-=[('°Sn)^S,(i/r^)]=^ (4.9) 

(for suitable C, c, functions of E), also summable in n: 
this bound can be taken from Appendix J in [2] where it 
is derived in detail in the text following Eq.(9.41). 
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